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ABSTRACT 

We report on the metal distribution in the intracluster medium around the radio galaxy 
4C+55.16 (z = 0.24) observed with the Chandra X-ray Observatory. The radial metallicity 
profile shows a dramatic change at 10 arcsec (~ 50 kpc) in radius from half solar to twice so- 
lar at inner radii. Also found was a plume-like feature located at ~ 3 arcsec to the south-west 
of the centre of the galaxy, which is mostly accounted for by a strong enhancement of iron L 
emission. The X-ray spectrum of the plume is characterized by the metal abundance pattern 
of Type la supernovae (SNela), i.e., large ratios of Fe to a-elements, with the iron metallicity 
being unusually high at 7.9^5^0 ^^^^^ (^^ P^^ ^^^^ error). How the plume has been formed 
is not entirely clear. The inhomogeneous iron distribution suggested in this cluster provides 
important clues to understanding the metal enrichment process of the cluster medium. 
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■ 1 INTRODUCTION 

, The metal abundance ratio pattern is useful for understanding the 

■ origin of heavy elements in the intracluster medium (ICM) because 
] of significantly different metal yields between Type la and II (plus 

lb and Ic) supernovae (SNela and SNell) (Renzini et al 1993). 
Overall, SNell appear to play a major role in enriching the ICM, 
as inferred from abundance measurements based on ASCA data 
(Mushotzky et al 1996; Fukazawa et al 1998). On the other hand, 
it has been suggested that an iron metallicity gradient is present in 
the ICM, especially of clusters with a cooling flow. Evidence for an 
increasing contribution from SNela at inner radii, inferred from an 
enahancement of iron abundance relative to a elements, has been 
reported for some clusters (Finoguenov, David & Ponman 2000; 
Fukazawa et al 2000; Kaastra et al 2001). However, the iron abun- 
dance generally remains significantly sub- solar and rarely exceeds 
1 solar even at the metallicity peaks in nearby clusters (e.g., De 
Grandi & Molendi 2001). 

Whereas the limited spatial resolution of the previous X-ray 
telescopes means that studies of a metallicity gradient in ICM are 
possible only for nearby clusters on a large scale, with the arc- 
second resolution of Chandra X-ray Obseratory (Weisskopf et al 
1996), we are now able to explore a cluster core at higher redshift in 
much more detail. In this Letter, we present spatially resolved high 
metallicity concentration in the central part of the cluster around 
4C+55.16. A comprehensive study of the cluster emission will be 
presented elsewhere. 

4C+55.16 is a compact powerful radio source residing in a 
large galaxy at a redshift of 0.240 (Pearson & Readhead 1981, 
1988; Whyborn et al 1985; Hutchings, Johnson & Pyke 1988). Re- 
cently, luminous cluster emission 10^^ erg s~^) around the radio 



galaxy has been recognized through ASCA and ROSAT HRI obser- 
vations (Iwasawa et al 1999). We assume Ho = 50 km s~^ Mpc~^ 
and qo = 0.5 throughout this paper. The angular scale is then 4.8 
kpc arcsec"^ for a redshift of 0.24. 



2 OBSERVATION AND DATA REDUCTION 

4C+55.16 was observed with the Chandra X-ray Observatory on 
2000 October 8. The galaxy was positioned on the ACIS-S3 detec- 
tor, which was operating in Faint mode at a focal plane temperature 
of —120° C. The background of the CCD chip did not show signifi- 
cant flare events but remained within 20 per cent fluctuation around 
the mean count rate most of the total 9.07 ks exposure. The data 
were reduced using software in CIAO 2.1.3 and the latest calibra- 
tion (caldb 2.7) released in August 2001. Despite improved quality 
of spectral fits relative to the previous calibration, the redshift in- 
ferred from the Fe K line is 0.2541° °^° , which we use in the spec- 
tral analysis but is higher than 0.240 from the optical spectroscopy 
of the cD galaxy (Lawrence et al 1996). 

Spectral analysis was performed using XSPEC version 11. 
Spectral parameters obtained from the analysis are presented with 
the 90 per cent confidence region for one parameter of interest 
throughout this paper, unless stated otherwise. 



3 RESULTS 

3.1 X-ray morphology of cluster core 

Adaptively smoothed 0.3-7 keV image of the cluster core obtained 
from the ACIS-S3 is shown in Fig. 1. An unresolved bright X-ray 
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Figure 1. Central part of the X-ray emission around 4C-I-55.16. The image 
size is 24 x 24 arcsec. The Chandra ACIS-S 0.3-7 keV band image has been 
adaptively smoothed with a 3cr gaussian kernel. Contours are drawn at 10 
equally-divided logarithmic intervals in counts per original pixel between 4 
and 40 per cent of the peak brightness occuring at the active nucleus of the 
radio galaxy. A plume-like feature is seen to the south-west. 

source is found at the position of the nucleus of the radio galaxy. 
The X-ray cluster emission is elongated along the major axis of 
the optical isophote of 4C-I-55.16 at a position angle of ^ 30°. A 
notable feature seen in the cluster core region is a plume-like X-ray 
excess extending by 5 arcsec to the SW of the nucleus. There 
is no obvious counterpart in the optical and near-infrared images 
(Hutchings et al 1988; C. Crawford, priv. comm.). 

3.2 Radial variation of the ICM properties 

The point source at the nucleus shows a hard X-ray spectrum, 
which can be attributed naturally to nonthermal emission from the 
active nucleus (Fig. 2). A power-law fit to the spectrum yields 
a photon index of 1.3 =b 0.3 with no significant excess absorp- 
tion above the Galactic value (A^h= 4.2 x 10^°cm~^, Dickey 
& Lockman 1990). The observed 0.5-2 keV and 2-7 keV fluxes 
of the nuclear source are 2.4 x 10~^^erg cm~^ s~^ and 5.7 x 
10~^^erg cm~^ s~^, respectively. The rest-frame 2-10 keV lumi- 
nosity is estimated to be 2 x 10^^ erg s~^. 

Also taken were spectra from eight annuli of radii of 1-2.5, 
2.5-5, 5-9, 9-15, 15-25, 25-40, 40-60, 60-100 arcsec. The back- 
ground data were taken from a region on the same detector where 
the cluster emission is negligible. Each spectrum was fitted by a 
single MEKAL model. For the metallicity measurement, we use 
solar abundance ratios taken from Anders & Grevesse (1989). The 
ACTS spectra and fitting results on temperature and metallicity are 
shown in Fig. 2 and Fig. 3. Excess absorption is not significant in 
most cases. However, it affects the temperature measurements and 
a few times 10^°cm~^ could still be within the calibration uncer- 
tainty. We therefore present results in Fig. 3 both for when Galactic 
absorption is assumed and variable excess absorption is allowed. 

The temperature of the cluster gas shows a monotonic de- 
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Figure 2. The ACIS-S spectra of the active nucleus and the ICM taken from 
annuli in radii of 1-2.5, 2.5-5, 5-9, 9-15, 15-25, 25-40, 40-60 and 60-100 
arcsec from top to bottom. Strong Fe L features are evident in the spectra 
from the inner three annuli. 



crease towards the centre. It is noted that flattening of the tempera- 
ture gradient is not seen up to 100 arcsec (480 kpc) in radius, which 
is rather unusual. A similarly high temperature is obtained for the 
outermost annulus when background is extracted from the blank 
sky data provided by Chandra Science Center. 

The metal abundance profile shows a dramatic increase at a 
radius around 10 arcsec. Within this radius, the metallicity signif- 
icantly exceeds solar whilst it is about half solar at the outer radii. 
A strong Fe L feature around 0.9 keV (in the observed frame), by 
which metallicity is primarily derived, is evident in the spectra from 
the inner three regions (see Fig. 2). Although individual spectra do 
not have sufficient counts in the Fe K band, the integrated spectrum 
within a radius of 9 arcsec (excluding the nuclear source) shows 
a clear strong Fe Ka line, supporting the high metallicity implied 
from the Fe L emission (Fig. 4). 
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Figure 3. The radial profiles of temperature and metallicity obtained from 
spectral fits to the eight annulus- spectra shown in Fig. 2 with a single 
MEKAL model. Results when the spectra are modified by Galactic absorp- 
tion (circles with solid-line error bars) and variable excess absorption is 
allowed (triangles with dotted-line error bars) are shown. The error bars are 
of Icr. The angular scale is 4.8 kpc arcsec"^ . 
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Figure 4. The ACIS-S spectrum for the central 9 arcsec region excluding 
the nuclear source, fitted with a single MEKAL model modified by Galactic 
absorption. A strong Fe K line is clearly detected at the observed energy 
around 5.4 keV, confirming the large metallicity inferred from the Fe L in 
the low energy band. The best-fit temperature and metallicity obtained from 
the fit are kT = 3.34tQ 3^ keV and 1.93tQ g2 ^^lar, respectively. 



3.3 Fe-L plume 

The X-ray plume located near the nucleus (see Fig. 1) is found to 
show particularly strong Fe L emission. The continuum- subtracted 
Fe L image, made from the 0.78-1.0 keV (0.95-1.25 keV in the 
galaxy frame) image, shows strong excess at the plume region 
(Fig. 5). Since the Fe L-free band (e.g., 1.8-5 keV) image shows 
a smooth elliptical isophoto centred on the galaxy nucleus with 
no particular excess at the plume, the plume feature is mostly ac- 
counted for by the excess Fe L emission. This also indicates that 
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Figure 5. The continuum- subtracted Fe-L emission image of the central 
part of the cluster. The 0.78-1.0 keV (in the observed frame) image has 
been adaptively binned so that each pixel contains at least 9 counts, and then 
the continuum for each pixel estimated from the 1.3-2 keV image was sub- 
tracted (the continuum level is 22 per cent at the Fe L peak). The colour bar 
shows surface brightness in counts per original pixel. Overlaid in contours 
is the adaptively smoothed 1.8-5 keV band image where the continuum is a 
major component and is free from Fe-L emission. The contours are drawn 
at 10 equally-divided logarithmic intervals of surface brightness from 2 per- 
cent of the peak brightness. Strong excess iron-L emission is seen at around 
~ 3 arcsec to the SW of the galaxy nucleus and accounts most of the plume 
feature, which is absent in the ICM represented by the continuum image. 
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Figure 6. The ACIS-S spectrum of the plume region fitted with the 
VMEKAL model with the SNIa abundance pattern (see Table 1). Note the 
strong Fe L emission around 0.9 keV. 



the ICM itself is not structured, unlike the disturbed X-ray mor- 
phology resulted from displaced ICM by radio jets seen in nearby 
clusters such as Hydra A (McNamara et al 2000) and Perseus clus- 
ter (Fabian et al 2000). In 4C-I-55.16, most of the radio emission 
originates from the VLB I- scale compact core and the small ra- 
dio lobe, aligned approximately perpendicular to the plume, carries 
only small power (Whyborn et al 1985). 



© 0000 RAS, MNRAS 000, 000-000 



4 K. Iwasawa et al 




Energy (keV) Z^^ (Solar) 



Figure 7. The left panel shows a ratio plot of the spectra from the C-shaped 
region (open squares) and the plume (filled circles) relative to the best-fit 
MEKAL model for the C-shaped region (Table 1), demonstrating excess Fe 
L emission of the plume spectrum. The plume data have been normalized 
so that the data excluding the Fe L band match the model. The data have 
been rebinned further for clarity. The right panel shows the curve for 
the metallicity of the plume spectrum when fitting with the SNIa abundance 
pattern (Table 1). This shows that the metallicity exceeds 1 solar at 3<t level. 



The spectrum of the Fe L plume (Fig. 6) is investigated. 
The spectral data were collected from a rectangular region of 
3 X 4.5arcsec with a position angle 0° , centred on the plume, which 
covers most of the SW quarter of the 1-5 arcsec annulus. We also 
take a spectrum from the rest of the annulus ("C-shaped region": 
approximately three-quarter of the 1-5 arcsec annulus covering po- 
sition angles from 280° to 180° to avoid the nucleus and plume), 
for comparison. The strong Fe L emission from the plume peaks at 
an energy of 0.89±0.02 keV, very similar to that (0.89±0.03 keV) 
of the C-shaped region spectrum, suggesting similar temperatures. 
Therefore the enhancement of Fe L emission seen in the plume is 
more consistent with a metallicity effect than the presence of cool 
gas. 

Fitting the C-shaped region spectrum by a single MEKAL 
model gives a good fit with a temperature of kT = 3.91q g keV 
and 3.412^ solar metallicity (Table 1). 

The metallicity inferred for the C-shaped region spectrum is 
already unusually high for ICM. The even stronger Fe L emission 
seen in the plume (Fig. 7) could be a result from selective enrich- 
ment of iron, presumably by SNela. Since our short exposure spec- 
trum of the plume (the integrated count is ^ 170) does not have suf- 
ficient quality for fitting emission lines from individual elements, 
we tried three different metal abundance patterns appropriate for 
a) Solar photosphere; b) SNII; and c) SNIa, in order to see which 
pattern reproduces the observed spectrum best. MEKAL is used for 
the solar abundance case, while the enrichment patterns by the two 
supernova types are implemented with the variable abundance ver- 
sion of MEKAL (VMEKAL), using the abundance ratios relative 
to the solar values given in Dupke & Arnaud (2001), based on the 
stellar yields for SNII (Nomoto et al 1997a) and for the deflagration 
model (W7) of SNIa (Nomoto, Thielmann & Yokoi 1984; Nomoto 
et al 1997b). 

The results of spectral fitting are shown in Table 1. The SNell 
type abundance pattern can be ruled out. The SNIa abundance pat- 
tern gives a significantly better fit over that of the solar abundance 
(Ax^ ~ 7 for 13 degrees of freedom for both fits). The SNIa abun- 
dance model gives a temperature kT 3.4 keV, similar to the 



Table 1. Results of spectral fits to the ACIS-S spectra of the plume 
and "C-shaped" regions. Three metal abundance patterns appropriate for 
(a) the solar photosphere; (b) SNIa; and (c) SNII were tested for the 
plume spectrum. All the fits include Galactic absorption. The number ra- 
tios of each element to iron, relative to the solar ratio, are as follows 
for the two types of supemovae, as given in Dupke & Arnaud (2001): 
(O,Mg)=0.035Fe, Ne=0.006Fe, (Si,S,Ar,Ca)=0.5Fe, Ni=4.8Fe for a SNIa; 
and (0,Mg,Si)=3.7Fe, (Ne,S)=2.5Fe, (Ar,Ca,Ni)=1.7Fe for a SNII. *K is 
the normalization of the model in unit of 10~^^ /{AttD^) J neUudV, 
where Dl is the luminosity distance of the source (4.7 x 10^^ cm), 
and riH are densities of electron and hydrogen, respectively in cm~^. The 
errors quoted are of 90 per cent confidence region for one parameter of in- 
terest. "'"This value is for the solar abundance pattern. Zpe = 2.3^ solar 
with a comparable is obtained when the SNIa pattern is employed. 
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surrounding C-shaped region, and 8 solar iron metallicity (as 
shown in Fig. 7, the iron metallicity exceeds 1 solar at 3(J confi- 
dence level). 

More complex models, e.g., two-temperature model and 
cooling-flow model, are not required by the data. The inclusion of 
extra components may be needed, which could alter the metallic- 
ity measurements if the data quality is much better, but it would 
usually push the metallicity value even higher. 

The model of the SNIa abundance pattern is also examined for 
the C-shaped region spectrum. It gives 2.3^ solar metallicity, but 
does not improve the fit over the solar abundance model, unlike for 
the plume spectrum (see Table 1). 



4 DISCUSSION 

We have found a striking increase in metallicity in the cluster core 
around 4C+55.16 through iron emission measurements. Metallic- 
ity increases of a lesser degree have been found earlier towards 
the centres of several nearby clusters (up to 0.5 solar, De Grandi 
& Molendi 2001). Centaurus and Virgo clusters show peak metal- 
hcity of - 1.5 solar (Ikebe et al 1999; Sanders & Fabian 2001) 
and ~ 0.8 solar (Bohringer et al 2001), respectively. Sersic 159- 
03 shows an iron-rich region (0.76 solar) displaced from the centre 
by ^ 30 kpc (Kaastra et al 2001). Our results for the higher red- 
shift cluster 4C+55.16 are more dramatic than any of these, with 
the cluster core being about twice solar abundance and the plume 
is apparently three times richer in iron. No mechanism other than 
a metallicity effect appears to be plausible: photoionization by the 
active nucleus, for instance, has a problem with orientation, i.e., the 
radio axis is perpendicular to the plume; and non-equilibrium ion- 
ization of the gas, which could be induced by a merger, would last 
for less than a million years and thus is unlikely to be observed in 
the absence of any optical signature of a merger. 

From the temperature, metallicity, and normalization of the 
VMEKAL thermal spectrum obtained from the fit to the plume 
spectrum with the SNIa abundance pattern, and its assumed vol- 
ume, V = 14.4 X 14.4 X 21.6 kpc^, we obtain gas densities. 
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rigas ~ 0.177 cm~^, and gas mass, Mgas — fJ^mpfigasV ~ 
1.2 X 1O^°M0, The iron mass, Mpe ^ 2.5 x lO^M©, the emis- 
sivity of the gas, e ^ 3.3 x 10~^^ ergs~^ cm~^, and its cooling 
time, tcooi = {3/2)kTngas/e p:^ 1.4 x 10^ yr are derived. 

The metal-enrichment of the plume region clearly needs iron- 
rich sources, such as SNela. With the iron yield per SNIa of 0.74 
M©, 3.4 X 10^ SNela are required to produce the iron in the plume. 
We discuss a few possiblities for the origin of the iron-rich plume. 

Frequent minor-merger events are expected in the cD galaxy. 
The plume could be the gas stripped from an infalling metal-rich 
galaxy, although iron needs to be selectively left behind. More 
likely is star formation taking place in the gas deposited by a 
merger. A short duration starburst results in a metallicity pattern 
of SNela after a significant number of SNela start to explode (e.g., 
Hamann & Ferland 1993; Yoshii, Tsujimoto & Kawara 1998). 
However, an immediate difficulty is the lack of excess in optical 
light there. The - 10^° M© or more of star formation for produc- 
ing the required number of SNela should be visible in blue light. 
There are two other potential problems, depending on evolutionary 
models of SNIa progenitors. Firstly, even if the progenitors were 
formed in a clump, they would have moved in the galactic potential 
and phase mixed within a timescale of order a few times the cross- 
ing time of the plume 5 x 10^ yr (its size, 14 kpc divided by the 
orbital velocity of say 300 km s~^). This is consistent only with 
the lifetime of currently disfavoured doubly-degenerate (DD) pro- 
genitors (Tutukov & Yungelson 1994) and single degenerate (SD) 
with a large mass (up to SM©) secondary star (Matteucci 1994), 
and not with that of SD binaries with a smaller mass secondary 
star (SNela only start to explode after 0.5 Gyr for an elliptical, 
Kobayashi et al 2000). Secondly, the time scale of SNela enrich- 
ment must match the short cooling time of the gas. For an elliptical 
galaxy with Lb ^ lO^^L© like 4C-h55.16, the peak SNIa rate, 
which lasts for less than 1 Gyr, is expected to be 0.1-1.5 SNela 
yr~^ (Kobayashi et al 2000; Matteucci 1994). With this range of 
SN rate, it takes 0.2-3 Gyr to produce 3.4 x 10^ SNela, hence the 
iron mass in the plume. Thus the cooling time of the plume requires 
a SN rate at the higher end of the expected range to be matched by 
the iron enrichment time scale. It should also be noted that a similar 
amount of iron is contained outside the plume within the central 25 
kpc. 

An alternative possibility is that the plume has been ejected 
from the centre of the galaxy. Gas enriched by SNela may be re- 
tained within the steep cluster potential until ejection takes place. 
However, in order for the gas not to cool and then drop out, some 
form of heating to combat the radiative cooling is required. SN 
heating can only partially offset the cooling; the thermal energy 
in the plume is {3/2)nekTV « 10^^ erg, which is equivalent to 
the total energy from the required 3.4 x 10^ SNela, assuming 10^^ 
erg per SN and no SN energy is radiated. Using the cooling time as 
a guide, the projected distance of 14 kpc between the plume peak 
and the centre of the galaxy requires the ejected gas clump to travel 
at a velocity of > 100 km s~^. Since no sign of a temperature 
rise in the ICM due to a shock is seen, the motion of the galaxy 
must be subsonic, i.e., the velocity should be less than the sound 
speed of ^ 600 km s~^ for the 3 keV intracluster gas. At a veloc- 
ity in a range of 100-600 km s~^ and the ICM density of ^ 0.1 
cm~^, whether ram pressure stripping works effectively is not triv- 
ial. Although a dynamical study of the cluster is required to be con- 
clusive, from an inspection of the X-ray and radio images, the cD 
galaxy appears to be lying at the centre of the cluster potential and 
the radio lobe which perpendicular to the plume does not show any 
sign of the radio-emitting plasma being swept back, indicating no 



evidence for fast motion of the cD galaxy (at least on the plane of 
the sky). One last possiblity is that the metal-rich plume is a fossil 
from when the cluster core was hotter. This requires that mixing of 
metals with surrounding gas is suppressed, much as conduction is 
reduced (e.g., Vikhlinin et al 2001; Ettori & Fabian 2000) 

In summary, the high iron metallicity found in the plume re- 
gion is unusual when compared with well-studied nearby clusters. 
Although the amount of iron could be supplied by SNela within the 
central galaxy, how it accumulated is not clear. 
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